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A thermal QCD Finite Energy Sum Rule (FESR) is used to obtain the 
temperature dependence of the axial-vector coupling of the nucleon, 
gA (T). We find that gA(T) is essentially independent of T, in the 
very wide range < T < 0.9 T c , where T c is the critical temperature. 
While gA at T = is (^-independent, it develops a q 2 dependence at 
finite temperature. We then obtain the mean square radius associated 
with gA and find that it diverges at T = T c , thus signalling quark de- 
confinement. As a byproduct, we study the temperature dependence 
of the Goldberger-Treiman relation. 



The possibility of creating a quark-gluon plasma in relativistic heavy ion 
collisions has sparked much interest in theoretical predictions for the on- 
set of this state |1[. In addition to the search for unambiguous processes 
signalling the formation of such a plasma, it is also important to under- 
stand the temperature behaviour of hadronic Green's functions and their 
associated parameters, viz. masses, widths, couplings, etc. The general con- 
sensus is that hadronic widths depend strongly on the temperature; in fact 
they are expected to diverge at some critical temperature T c , thus signalling 
quark-gluon deconfmement || (hadronic widths are to be understood, in this 
context, as absorption coefficients determined by the imaginary parts of two- 
point functions). Thermal three-point functions also provide independent 
evidence for this phase transition, as the mean square radii happen to in- 
crease with increasing temperature, becoming infinite at T = T c ||. 
A recent investigation of the thermal behaviour of the pion-nucleon cou- 
pling, in the framework of both the linear sigma model and QCD sum rules 
0, showed that as the temperature approaches T c , g n NN(T) vanishes, while 
the associated radius diverges. Both g^NNiT) and (r% NN )(T) may thus be 
interpreted as signals for the deconfinement phase transition. In this work 
we shall determine the temperature behaviour of the axial-vector coupling 
constant of the nucleon g& = <m(<? 2 = 0), and the associated radius, using 
the method of thermal QCD sum rules ||; specifically, the leading dimen- 
sion Finite Energy Sum Rule (FESR). However, we shall first discuss our own 
determination of g^ at T = 0, as previous QCD sum rule determinations, 
dating back many years 0, were the subject of some controversy. We find 
it possible to reproduce the experimental value of Qa at T = 0, which then 
serves to normalize the finite temperature results. Finally, as a byproduct, 
we shall use this result to determine the behaviour at finite temperature of 
the SU{2) L x SU(2) R Goldberger-Treiman relation 

U( t )9ttNn(T) _ 
M N (T)g A (T) i + A ^ • W 

In this relation, f n (T) is known up to T = T c J7|, where it vanishes, g n NN(T) 
behaves qualitatively similarly |4j], and M^{T) is essentially constant up to 
T = T c ||-0. The question is then, how big are the thermal corrections to 
this relation, A 7r (T) (normalized to A 7r (0) = 0). An equally important chiral- 
symmetry relation, the Gell-Mann, Oakes and Renner relation (GMOR), has 
recently been investigated in the framework of thermal chiral perturbation 
theory [HJ and QCD sum rules |ll|. There is excellent numerical agreement 



between both results, indicating that temperature corrections to the GMOR 
relation are rather small. It should be kept in mind that a comparison be- 
tween thermal QCD sum rules results and those from effective theories at 
finite temperature, e.g. sigma model, chiral perturbation theory, etc., must 
necessarily be done numerically. The fields involved in the former technique 
are those of the quarks and gluons, while those of the latter framework are 
purely hadronic. As a result, expansions in powers of the temperature do 
not necessarily need to match order by order because the coefficients in these 
expansions will involve different types of parameters. However, numerical re- 
sults from both techniques should agree, at least within the range of validity 
of the low temperature expansion in effective theories (QCD sum rules are 
in principle valid across the whole range of temperatures). This is precisely 
what happens with the two analyses of the GMOR relation mentioned above. 

We begin by considering the three-point function 

n M (p,j/,g) = t* [J d t xdfy(0\T(r h (x)Mv)l«W)\Q)^'~™ ) , (2) 

where the charged axial vector current is given by: A^x) = u{x)^^d{x), 
while the interpolating currents of the proton and neutron are chosen as [T2j 

Vp(x) = e abc [u a {x)C^^u b {x)] Ylbd c , 

Vn(x) = -e ahc [d a (x)C lfl d b (x)] 7 M 75« C . (3) 

The axial- vector coupling of the nucleon, 5u(? 2 ), is defined through 

(N(p 2 )\A,(0)\N( Pl )) = u(p 2 ) [^J 5 g A (q 2 ) + g M 75 h A {q 2 )] u( Pl ) , (4) 

with q^ = (j>2 — Pi)fi- The coupling of the interpolating currents, Eq.(3), to 
the nucleon is 

(0\ V (0)\N(p)) = \ N u(p) . (5) 

Inserting a complete set of intermediate nucleon states into Eq.(2), one ob- 
tains the hadronic representation 

n„(p,p', q) = ( p 2 _ M 2^,2 _ M 2 n) (/ + M N )TM+ M N ) , (6) 

where 

Tf, = [7^75 gA(q 2 ) + q^lb h A (q 2 )] , (7) 
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and the following expansion holds 

(/ + M N )T^+ M N ) = g A (q 2 ) [-2 p"p^ 5 + / ^75 + (j+ f)l^M N 

-2 Pli -y 5 M N + M 2 Nlfll5 ] 
+ h A (q 2 ) [- p" ^+ 4M N + A4] g M 75. (8) 

Since we are only interested in qa, we need to extract tensor structures which 
are not multiplied by Ha] a suitable candidate being the structure 

(/+ J*)7m75 • (9) 

The relevant term of the imaginary part of the (hadronic) correlator is then 

Im iys,s',g 2 )| H AD = -\ 2 N g A (q 2 )M N 7r 2 5(s - M 2 N )5(s' - Af£)(/+ ^) 7/ ,75 

+ e(s - s )Q(s' - s' ) Im n M (s, s', g 2 )| QCD , (10) 

where s = p 2 , s' = p' 2 , and we have added the hadronic continuum, modelled 
by perturbative QCD, starting at thresholds s = so and s' = s' . Considering 
the contribution to the correlator from perturbative QCD, we obtain 

n M (p,p',g)|p Q cD = -24i 2 / X 2 -f a S F (k 2 )-f P S F (k 1 - q)^^ 

xSpikJ-fSptf - h - fchV (11) 

Taking the imaginary part of this expression, and evaluating the integrals, 
it turns out that there are no terms proportional to the tensor structure of 
Eq.(9). Turning to the non-perturbative part, we find the quark condensate 
contribution to the correlator to be 



n>,p',g)| Q cD = 2t 3 {dd) 



d 4 k 

■laS F (k)7[W»S F (q)7 a S F (p' -k- g)7 /3 7s 



+ (uu) 



(2vr) 4 
f -0j- 4 l a lpS F (k - q) lfl S F (kh a S F (p> - k)^ l5 

f-0^ 4 l a S F (p' - k) lp S F (k - qh,S F (kh a l\ 

r -0^1*S F (k) lp S F (q) 1 ,rS F (p' - k)^ l5 (12) 



Taking the imaginary part, and keeping only terms proportional to the rele- 
vant tensor structure Eq.(9), and which are non- vanishing in the limit q 2 — > 0, 
we obtain, after assuming (uu) ~ (dd) = 



Im n M (p,p',g)| Q cD = Y2^(<?9) (/+ /)7mT5- (13) 

Next, using Cauchy's theorem, and assuming quark-hadron duality, the low- 
est dimensional FESR for qa reads 

rs rs' /-so i>s' 

/ / ds ds' ImTl fM (s,s')HAD= / / ds ds Im II M (s, s')q C d • (14) 
Jo Jo Jo Jo 

From this FESR one then obtains the relation 

Sos' (qq) . , 

* A = ~12^Ap^' (15) 

At first sight, this result hardly looks like a prediction for gA, since so, s' , 
and Aat are a-priori unknown. However, since the double dispersion in p 2 = s 
and p' 2 = s', used in obtaining Eq.(15), refers to the nucleonic legs of the 
three-point function, it is reasonable to set Sq = s' . At the same time, a 
QCD FESR analysis of the two-point function involving nucleonic currents 
@] yields the following relations 

A^ = -A_ , \* N M N = -M s 2 (16) 

N 192tt 4 ' N 8tt 2 ° ' v ; 

where, in principle, the numerical value of the asymptotic freedom threshold 
So does not have to be the same as that in Eq.(15). In fact, if one were to 
assume them to be equal, then Eqs.(13) and (14) would imply qa = 8/127T, a 
value far too small. Without any attempt at extracting a precision value of 
gA-, it is rewarding, though, to find that the experimental value g^ = 1-26 can 
be reproduced in this framework if A^ ~ 3.1 x 1CT 4 GeV 6 , s = 3.7 GeV 2 , 
and the standard value (qq) = —0.01, are used in Eq. (15). These values of 
X 2 N and So are close enough to those resulting from the two-point function 
channel, at least for the purpose of the present work, which is to obtain the 
temperature dependence of gA and its mean square radius, rather than a 
prediction for gA (T = 0). 



The finite temperature corrections to g A are obtained by inserting the 
thermal Dolan and Jackiw |13|] propagators, and allowing for the temperature 



variation of (qq), Xn and So- For (qq)r and Xn(T) we shall use the results 
of and of ||, respectively. The temperature dependence of sq was first 



obtained in ||14| , and later improved in |p^| . It turns out that for a wide 
range of temperatures not too close to T c , say T < 0.8T C , the following 
scaling relation holds to a good approximation 

The appropriate contribution to the thermally corrected QCD spectral func- 
tion becomes 

Im U,(p,p',q) = M(^+ / )7m75 [ffaT) + /(p ' 5T) ] f ( i 8) 

where 

f( P , T) = £ ta [l - n F (^m ) - „, (!» + B*l)] , (19) 

with np(x) = (1 + e x ) _1 , and f(p',T) is similarly defined. Finally, we obtain 
the sum rule for g^ at finite temperature: 

(w) i r o(T) r° (T) 

9A ^ = -^WJT ds ds'if(p,T) + f( P ',T)]. (20) 

487T° A N M N J J 

In order to evaluate the integrals one needs to choose a specific frame, for 
example the (rest) frame p = 0. In this case, the components of the four 
vectors p and p' may be expressed in terms of s, s' and q 2 . Other choices 
of frames give essentially the same results. A numerical evaluation of gA{T) 
is presented in Fig. 1. As can be seen from this figure, g A is basically T- 
independent, and it clearly does not vanish as the critical temperature is 
approached. In this sense, gA does not represent a signal for the decon- 
finement phase transition. We turn now to the mean square radius (t a )t 
associated with g A , and defined as 

(r 2 A ) T = Q-^lng A (q 2 ,T)\ q 2 =0 . (21) 



This radius is non-zero at finite temperature due to the g 2 -dependence of the 
arguments of the thermal Fermi factors. After evaluating the logarithmic 
derivative of Eq.(20) one obtains 

f /-so rs'o ) -1 rso(T) ps' (T) rl 

(r\)T = { / dsds'[f(p,T) + f(p',T)]\ / ds ds' dx 
Jo Jo Jo Jo J-l 



6 / p' \ p' +\p'\* 



Up 



' Vo + W '\ x 



1 2 



(22) 



This is plotted in fig. 2, which shows that the radius diverges as the critical 
temperature is approached. This kind of behaviour has been obtained previ- 
ously for other radii ||, ||, and it may be interpreted as (analytic) evidence 
for quark deconfinement. 

Finally, we can use our result for gA at non-zero temperature to evaluate 
the validity of the GTR, Eq.(l). Results for the mass of the nucleon show 
that it has very little variation with temperature, and so we shall assume 
that it is constant ||-||]. Using the result of |7| for f n at finite temperature, 
together with our previous results for g V NN(T) ||], and our current result for 
gA{T), we can determine the thermal correction to the GTR, A 7r (T) defined 
in Eq.(l). In fig. 3 we present a plot of 1 + A^T) against T/T c , which indi- 
cates that the GTR is approximately correct until about T ~ 0.9 T c , where 
it breaks down. 
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Figure Captions 

Figure 1. The coupling qa(T), from Eq. (20), as a function of T/T c . 

Figure 2. The temperature dependence of the mean square radius, Eq.(22). 

Figure 3. Deviation from the Goldberger-Treiman relation, Eq.(l), as a func- 
tion of of T/T c . 
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